The semi-analytic planetary theory VSOP (French: Variations Séculaires des Orbites Planétaires) is a mathematical model describing long-term changes (secular variation) in the orbits of the planets Mercury to Neptune. The earliest modern scientific model considered only the gravitational attraction between the Sun and each planet, with the resulting orbits being unvarying Keplerian ellipses. In reality, all the planets exert slight forces on each other, causing slow changes in the shape and orientation of these ellipses. Increasingly complex analytical models have been made of these deviations, as well as efficient and accurate numerical approximation methods.

VSOP was developed and is maintained (updated with the latest data) by the scientists at the Bureau des Longitudes in Paris. The first version, VSOP82, computed only the orbital elements at any moment. An updated version, VSOP87, computed the positions of the planets directly at any moment, as well as their orbital elements with improved accuracy.